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We study the top quark portal dominated dark matter interactions, and its implications for the 
gamma ray line searches. In this picture, the dark matter interactions with photons and gluons 
are loop induced by the axial anomaly of the top quark current. We show there can be a natural 
suppression of the tree-level annihilation of dark matter, and the photon channel in turn has a 
substantial rate when the main annihilation proceeds into gluons. We observe a competition between 
the indirect detection of gamma ray line and the search with monojet plus missing energy events at 
LHC, and the 7 TeV data already set an upper bound of ~ 10~ 28 cm 3 s _1 on the photonic annihilation 
cross section. This upper limit is compatible with a thermal WIMP scenario. 



Introduction. The existence of dark matter (DM) in 
the universe has been established by various cosmologi- 
cal observations, yet its identity has not been uncovered. 
Many experiments are now going on in order to directly 
or indirectly detect its trace. The indirect detection us- 
ing the monochromatic cosmic gamma ray could serve 
as a clear evidence for DM being from particle physics 
origin. Since the DM is electric neutral, its annihilation 
into two photons must happen at loop level, naively sup- 
pressed by a factor of (a/ir) 2 compared with the tree level 
annihilation to charged final states. For thermal weakly 
interacting massive particle (WIMP) picture whose to- 
tal annihilation cross section is around a picobarn, the 
induced two photon annihilation rate is far below the 
current Fermi LAT sensitivity. 

Early this year, there were several analysis [1-3] of the 
4-year Fermi data showing the positive signatures in the 
gamma line search from the galactic center, with energy 
130 GeV and a cross section 1.3 x 10~ 27 cm 3 s~ 1 if origi- 
nating from DM annihilation [4]. This is only one order 
of magnitude below the thermal cross section. A more re- 
cent data reprocessing performed by the Fermi-LAT col- 
laboration shows a similar line-like feature which slightly 
shifts to 135 GeV. While the significance of this excess 
remains at ~ 3 er, it seems at present one cannot exclude 
the uncertainties from systematics or the earth limb pho- 
ton background [5] . On the conservative side, the global 
analysis of Fermi data offers an upper bound on the DM 
to two photon annihilation cross section, which is around 
10~ 27 cm 3 s -1 for photon energy 130 GeV [6]. 

The above hint has inspired a plethora of theoretical 
studies [7, 8, 14] on DM models that can give an "en- 
hanced" gamma ray feature from annihilation. Questions 
need to be addressed include: whether this signal can 
be reconciled with thermal DM paradigm; what are the 
correlated phenomena and in turn their constraints [9- 
11]. Possible answers to the former question involve the 
suppression or elimination of direct annihilation to light 
charged standard model (SM) fermions. In many cases, 
additional states other than the DM itself are introduced, 
playing the role of co-annihilator or intermediate states 
(real or virtual) in the annihilation. 



In this paper, we investigate a case when DM \ couples 
predominantly with the top quark [12-14], and assume 
the couplings to other SM fermions or bosons are neg- 
ligibly small throughout the discussion. For the bulk of 
this work, we consider the case when the DM is lighter 
than the top quark. In this case, it cannot annihilate 
into both on-shell top quark and antiquark, and we no- 
tice a fact that the three-body annihilation threshold 
mt + mw + mb ~ 256 GeV, is only a few GeV below twice 
of the suspected gamma line energy in Fermi (~ 130 - 
135 GeV), or twice of the corresponding DM mass. This 
indicates additional final state phase space suppression 
for the process XX ~^ tW~b or tW + b. Meanwhile, the 
DM can also annihilate into two photons via the top 
quark loop. The ratio between the two cross sections 
can be estimated as 



crv xx ^.twb J 8m? ' 

where the factor Air stands for a generic ratio of two- and 
three-body final state phase spaces, and 8m ~ 2m x —m t — 
mw — mb provides additional suppression when the later 
is being kinematically squeezed. For m x ~ 130 GeV and 
8m ~ few GeV, the ratio can be as large as 0.1. It serves 
as one motivation to study possible enhanced gamma ray 
line signals with this setup. Another advantage of the top 
quark portal is that the continuum gamma ray constraint 
turns out to be less severe [11]. 

A particular feature of the top quark loop induced an- 
nihilations into 77 and gg discussed here is that, they al- 
ways proceed through the axial anomaly of the top quark 
current. For the mass range m x < 130 GeV, the XX 99 
channel will dominate over the three-body one and lead 
to the largest ratio of photonic to the total annihilation 
cross section, which is ~0.4%. The DM mass indicated 
from the Fermi data lies within this window. 

The main consequences of the above picture are as 
follows. Such top quark portal being behind the Fermi 
gamma ray line implies a sizable effective coupling of DM 
to gluons. The total annihilation cross section exceeds 
that required for the thermal relic density. This case is 
severely constrained by LHC measurement of missing en- 
ergy plus monojet events and has been excluded using the 



7TeV data. On the other hand, the thermal WIMP case 
is still allowed by LHC data. The two photon annihila- 
tion cross section is about one order of magnitude smaller 
than that needed to explain Fermi, but large enough to 
be probed with more data and future experimental sen- 
sitivity [15]. We foresee an intimate interplay between 
DM collider searches and indirect detections. 

Top Quark Portal and DM Annihilations. In 

this work, we consider the DM as a Dirac fermion and 
write down its effective couplings with the top quark. In 
order to avoid large coupling to the lighter quarks, such 
as the bottom, and preserve the SM gauge symmetry, we 
choose only the coupling to the right-handed top quark. 
The viable operator at lowest dimensions is [16] 

Mortal = ^ frV (° + 7ti)x] [«7"(1 + 75)*] , (2) 

where A is the cutoff of the effective interaction and a is 
a real number. 

As sketched in the introduction, the annihilation of 
DM has the following channels (see also Fig. 1 and 3) 

( tW~b, iW+b 
XX -> I 77, lZ (3) 
[ 99 ■ 

The annihilation to 77, 7Z and gg all happen via virtual 
top quark loop, which we calculate in details below. 




FIG. 1. DM annihilation channels. Left: tree-level three- 
body process, kinematically suppressed in the phase space. 
Right: top quark loop induced annihilation into two photons. 

Annihilation XX ~^ tWb. 

For the three-body annihilation, we obtain the amplitude 
square 

x {{a + l) 2 (p 2 ■ p 3 ) [M^(jn ■ P5 ) + 2(pi • Pi){ P i ■ P5 )] 
+ (a - lf{ Pl • p 3 ) [M^(p 2 ■ p 5 ) + 2(p 2 ■ p 4 )(p 4 • p 5 )] 
+ (a 2 - l)m 2 [M^(ps • p 5 ) + 2(p 3 • Pi)(Pi ■ p 6 )] } , (4) 

where t, W b are labeled as 3, 4, 5 respectively, the pre- 
factor 2 takes into account of the charge-conjugation pro- 
cess. Here s « 4(m x + T/) 2 or 4m 2 during the DM freeze 
out and today, respectively. In the numerical calculation, 
we set the freeze out temperature Tf ~ m x /25, in order 
to approximate the thermal average. 

We integrate over the general three-body final state 
phase space as given in the appendix [17]. The resulting 
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FIG. 2. Cross section for the three-body annihilation \X ™* 
tWb during freeze-out (Tf ~ m x /25, which corresponds to ve- 
locity v f m 0.3), in units of cm 3 /s. Here we set a = 0. Notice 
such annihilation will not happen today if m x < 128 GeV. 

cross section during the freeze out is shown in Fig. 2. As 
one can see, the three-body annihilation can be highly 
suppressed by kinematics in the final state phase space, 
when the DM mass approaches the threshold for the 
t, W, b final states from above. The annihilation of DM 
today is more suppressed or even forbidden due to much 
lower velocity. 

In view of this suppression, it is meaningful to examine 
the other loop level processes. 

Annihilation XX ~^ 77- 

We first examine the DM annihilation into a pair of pho- 
tons, which leads to a spectacular line on top of the 
gamma ray spectrum. The DM annihilation into pho- 
tons can take place when the top quark lines are closed 
into a loop (see Fig. 1). When the tree- level annihilation 
is suppressed as shown above, the loop-level annihilation 
into photons could be more notable. 

Using the interaction Eq. (2), the annihilation of DM 
into two photons arises from the axial current £7^75*, 
evaluated inside the triangle loop in Fig. 1. The cor- 
responding contribution from a vector current vanishes 
due to the charge-conjugation invariance. Following the 
generic calculation of axial anomaly [18, 19], the finite 
amplitude is 

■*x*-rr = ^^^3(fci,fe)f; x (p 2 )75%(pi) £ ' £lfc2eie2 , 

(5) 

where on-shcll photon conditions have been implemented, 
&i,2 are the photon momenta, and the factor 3 stands 
for the color factor. If the top quark is integrated 
out, the effective operator for the annihilation becomes 
m x x75%FF. The annihilation of DM to photons hap- 
pens only through the axial current anomaly, and only 
the axial part of the dark matter current contributes. 
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This also agrees with the argument [20] that, the process 
vv — ¥ 77 vanishes if the initial neutrinos are massless. 

Working in the center-of-mass frame, the form factor 
^4-3 (fei, fe) simplifies to 
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The resulting cross section is 
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As a sample value, for m x = 130 GeV and A = 85 GeV, 
we get erv xX _ i . 77 = 1.21 x 10 -2 ' cm 3 S -1 , which is close 
to the central value that explains the Fermi line excess. 
Notice a fairly low cut-off scale A is required here, which 
indicates the higher scale physics may already be accessi- 
ble. One possible UV completion for the effective portal 
could be the t-channel exchange of a colored and Zi odd 
particle, which based on the structure of Eq. (2) and 
upon the Fierz transformation, corresponds to a = 1 for 
a vector particle, and a = — 1 for a scalar one. The di- 
rect pair production and decay of such particle at LHC 
will result in the excess of it + $t events. However, if 
the mass of such particle is very close to the sum of top 
quark and DM masses, it becomes more difficult to peel 
from the SM background [21]. 




FIG. 3. DM annihilate into gluons. 
Annihilation XX ~^ 99- 

The same top quark loop also induces the annihilation of 
XX ~^ 99 1 shown in Fig. 3, which is a potentially dominat- 
ing channel during the DM freeze out. Its cross section 
can be rescaled from the two photon one, even without 
integrating out the top quark. With the electric charge 
Qt = 2/3, the ratio is fixed 
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This implies an upper bound on the ratio 
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(8) 



(9) 



which is saturated when the DM mass is such that the 
annihilation into gluons dominates over the three-body 
one, av tot w av xx ^ gg . 

In Fig. 4, we compare the two potentially dominant an- 
nihilation channels, XX ~^ tWb and XX ~~ ^ 99- To the left 
of the thick curve, the DM annihilation during freeze out 
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FIG. 4. Ratio of DM annihilation cross sections during 
the freeze out epoch, av x ^^ t wt/o'v x ^^^ 1 . Since the ratio 
o-« xx _ >ss /o-« xx _ >77 = 258 is fixed, we find the DM annihila- 
tion is dominated by the (anomaly loop) induced annihilation 
into gluons \X ~ * 99 m the region to the left of the thick curve. 



and today is controlled by the loop (axial anomaly) in- 
duced processes. For sufficiently small |et|, this regime of 
parameter space covers the ~ 130 GeV DM mass range as 
hinted by Fermi. In this case, our scenario could provide 
a minimal realization of the picture in [8], with a rela- 
tively smaller electric charge of the top quark though, 
and without introducing exotic heavy quarks. 

On the other hand, if the dark matter is a thermal 
WIMP, from the upper bound Eq. (9) , the maximal pho- 
ton annihilation is on the order of ~ 10~ 28 cm 3 s -1 . 

Annihilation XX ~^ 

It has been argued [22] that the two photon annihila- 
tion of DM could be accompanied with the annihilation 
to Z7, which is also the case here. This leads to a sec- 
ond line in the gamma ray spectrum at 114 (120) GeV 
for m x — 130 (135) GeV. Here XX ~^ %1 is generated 
from the same top quark loop, but the axial current can 
be placed at either the Z vertex or the DM portal. We 
find the following ratio among different annihilation cross 
sections 



avz-y ■ avyj w 0.76 : 1 



(10) 



corresponding to m x — 130 GeV and a = 0. The subse- 
quent decays of Z-boson into charged fermions near the 
center of galaxy induces synchrotron radiation, which is 
within the sensitivity of the galactic radio telescopes [10]. 

Collider Implications. As shown in Eq. (8), there 
is a sizable ratio between the DM annihilation cross sec- 
tion, or the effective coupling to gluons compared to pho- 
tons. This means if the gamma ray line predicted in this 
scenario is large enough to be probed presently, the cor- 
responding gluon coupling will be constrained by other 
experiments, e.g., from the anti-proton flux in the cosmic 
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ray to hadron colliders. While the former is often sub- 
jected to larger uncertainties among the DM halo profiles 
and the cosmic propagation models, the search using the 
monojet plus missing energy events at LHC could give 
a relevant and more solid constraint, which we focus on 
below. 

In this model, the main interactions of DM relevant 
for the collider study are those with gluons. Similar 
to Eq. (5), assuming the form factor could be approx- 
imated by a contact one in the heavy top quark limit, 
^3(^1, k 2 ) — > 2-n- 2 / (3rn 2 ) , the effective operator becomes 
X15XGG. The parton level process is gg — > xX+J e ^ with 
the additional jet emitted directly from the effective ver- 
tex, i.e., the top quark loop, or from initial states. 

We extract the current LHC constraints on the top 
quark portal based on the following considerations. The 
analysis in [23] shows that the two effective operators 
XxGG and XlixGG share similar monojet constraints. 
The operator a s /(4Alf )xxGG has already been con- 
strained by the 7 TeV LHC data [24], which is M* > 
300 GeV for DM mass around 130 GeV. In the contact 
interaction approximation, the corresponding Wilson co- 
efficient is 



a s m x 



4M| 37rA 2 m t 2 ' 
and the ATLAS limit gets translated into 
A > 200 GeV . 



(11) 



(12) 



A more precise calculation of parton level cross section 
keeping the explicit form factor A^(k\,k2) will slightly 
enhance the predicted signal, and leads to a relatively 
stronger bound on A. 

Both the DM annihilations to gg, 77 and the mono- 
jet constraint depend on a single scale A. Therefore the 
resulting collider constraint is closely connected to the 
indirect detection with gamma ray lines. From the AT- 
LAS result, we obtain the following implications for the 
top quark mediated DM scenario, which serve as the main 
message of this work. 

• The lower bound on A from monojet search has 
already exclude the parameter space that could ac- 
commodate the gamma ray line excess at Fermi. 
The is mainly due to the fixed and large hierar- 
chy between the gg, 77 effective couplings to DM 
induced by the top quark loop. 

• In the window where \X ~^ 99 dominates the to- 
tal annihilation rate, the case of DM being a ther- 
mal candidate is still but marginally consistent with 
the current monojet constraint. The corresponding 
XX ~^ 77 cross section is < 10 _28 cm 3 s _1 . It will 
keep being tested by the upcoming data from both 
collider and indirect detection experiments [15]. If 
LHC soon sees excess in the monojet events from 



DM, then a gamma ray line signature should also 
be around the corner, and vice versa. There is an 
interesting competition. 

For DM much heavier than 130 GeV, and the monojet 
constraint gets weakened. Moreover, when the tree-level 
annihilation dominates the freeze out, one needs to take 
higher value of A for it to be thermal. In this case, the 
contact top quark portal may be probed by other chan- 
nels such as ti + where a pair of DM is produced in 
together with top-anti-top pair [12]. 

Direct Detection. In this section, we briefly sketch 
the direct detection of DM from the above top quark 
portal. At the energy scale of DM-nucleon scattering, 
the top quark is heavy and can be integrated out. 

First, the anomaly loop induces a set of non- 
conventional effective gluonic operators for the direct de- 
tections. The effective interactions are, 



£eff — 
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m x XH5XG^G^ +xri,x(D a G a ")G p 



(13) 



where G^ v = \t ilvp ' J G pa . In the second operator, the 
derivative on G ap has been promoted to a covariant one 
in respect of gauge invariance. It vanishes when con- 
sidering DM annihilating into gluons due to the on-shell 
condition, but will contribute here to the matrix element 
between nucleons. Both operators lead to spin-dependent 
scatterings. 

The nucleon matrix element for the first gluonic oper- 
ator is 



32tt 



(N\G pu G» v \N) = ^- g ^m N u N i 75 u N , (14) 

where g^ « 0.36 [25], and N ~ p,n stands for pro- 
ton and neutron, respectively. The resulting amplitude 
(^x75 M x)(^Jv75 u ^)' an< l m t urn t ne differential cross sec- 
tion dtj/dEn, is suppressed by the momentum transfer of 
the DM-nucleon scattering, which is apparent after the 
non-relativistic expansion. For the choice of parameter 
A that can account for the monochromatic photon cross 
section within the Fermi LAT sensitivity, the correspond- 
ing spin-dependent cross section is found to be less than 
10~ 20 pb, far below the recent direct detection [27, 28] 
and neutrino telescope [29] bounds. 

The second operator in Eq. (13) is equivalent to, by 
using the equation of motion, 

(N\(D a Cr p )G PM \N) = J2( N \93QfGp P Yq f \N) 
f 



— ^ 2f2f m 'N S fj, 



(15) 



It is a spin-dependent twist-4 operator, whose ma- 
trix element is proportional to the nucleon polarization 
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Sfj, [26]. As a rough estimate of the coefficients f£j, as- 
suming up/down quark flavor dominance inside a pro- 
ton/neutron state, one can infer the sum of coefficients 
from the one (/| = J2f e qfl°f) quoted in [26] and esti- 
mate 



fP , fP , f 

J2u ' J 2d ' J'. 



(9/4)/| 



0.1 



(16) 



Therefore, such operator also leads to a spin-dependent 
DM-nucleon interaction, with effective operators at nu- 
cleon level, 2£/v(x7' i 75x)(-^ s A»-^0- The coefficient is 

The resulting spin-dependent direct detection cross sec- 
tion on a proton is 



XP 
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Vf 2 



(18) 



where zx is the reduced mass \i = m p m x /(m p + m x ). 
Using the same values of parameters as above, the corre- 
sponding spin-dependent cross section on a proton is less 
than 10 _10 pb, also well below the current experimental 
limits. 

The reason for this suppression can also be understood. 
This effective operator is dimension eight, suppressed by 
l/(A 2 m 2 ). The nucleon matrix element of the gluonic 
part of operator is typically of order m N . This leads to 
an additional suppression of (mjv/A) in the final cross 
section, compared to the usual WIMP case. A possible 
way to enhance such anomalous top quark current in- 
duced DM scattering to probeable level could be trading 
the effective scale suppression 1/A 2 for sub-GeV force 
mediators. However, this would require further model 
building and the existence of new light states could in 
turn modify the above DM freeze out processes and de- 
tections. We leave a systematic study of it for a future 
work. 

It turns out the most promising direct detection of such 
DM from the top quark portal is the exchange of a sin- 
gle photon in the t-channel [30], which is spin indepen- 
dent. Although this operator is suppressed by the small 
electromagnetic coupling, the photon exchange yields a 
dimension 6 operator, therefore relatively enhanced by a 
larger nucleon matrix element, compared to the gluonic 
operators. 

Without specifying the UV complete theory, we could 
estimate the corresponding cross section averaged per nu- 
cleon as 



will continue to be tested with more data, and the future 
Xenon IT [32] and LUX [33] experiments. 

Conclusion and Outlook. To summarize, we have 
studied a scenario where the DM interacts with the SM 
sector predominantly via the top quark portal. This pic- 
ture is partly motivated by the recent hint of gamma 
ray line excess in the Fermi data. The tree level anni- 
hilation with one top quark off-shell can be kinemati- 
cally suppressed, for the DM mass near 130 GeV. There 
is a window where the DM annihilation is dominated by 
axial current anomaly induced coupling to gluons (pho- 
tons). In this window, the annihilation into photons is 
at a percent level of the total cross section. We show the 
current 7TeV LHC monojet data is already in tension 
with the gamma ray hint at Fermi in this picture. In 
contrast, a thermal WIMP is marginally allowed by the 
LHC data, and the corresponding monochromatic pho- 
ton flux is within the sensitivity of current and planned 
indirect detection experiments. 

It will be interesting to see the future interplay between 
the DM searches at colliders and indirect detections. 
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Appendix: three-body annihilation. Here we 
give the general final-state phase space integration for 
the three-body annihilation cross section %(l)x(2) — > 
t(3)W(4)6(5), following the notation in [17], 



(p§)= 



dpi / dp\ 



dcost 



d I" 4 '" 

n ^8(2^)4$ 



(20) 



where 6, r] is the angle between p 5 and the z, y 
axis, respectively, and $ = 4^/1 |ki|. All the the 
relevant momentum products in the amplitude square 
can be expressed as functions of independent vari- 
ables p®, p®, 9, rj. The integral limits are (p^)z 
- [(m t + M w f - m 2 ]/(2^s), and (p°) 

j_ 

2r 



max 
.O^max 
min 



with a 



o(t + m + m_) ± y/ (p^) 2 — "m\\J ' (j ~ ™+)(t — m 2 _) 
fs—p®, t = a 2 — (pi) 2 +m 2 , m± = m t ±M n 



» 2 Z 2 (3aQ t Y 
it A 2 \2irA 2 J 



(19) 



in the case of a = —1 [30]. The cross section is a few 
10~ 45 cm 2 for the xenon target and A ~ 200 GeV (fol- 
lowing the LHC monojet constraint). This is close to 
the latest Xenon 100 limit [31], and the parameter space 
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